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Abstract 

Using the formalism of Bethe-Salpeter equation (BSE) the Higgsonium bound state is studied. 
The condition for formation of Higgsonium bound states are discussed in SM and in the simple 
extension of. 
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I. INTRODUCTION 



The Standard Model (SM) of particle physics involves very minimalistic idea of elec- 
troweak symmetry breaking scenario. Consequently it results with the only one real scalar 
field and in fact its last directly unobserved particle- the Higgs boson. However, being 
inspired of family pattern of SM fermionic sector, it is quite natural to consider the ex- 
tensions of SM with more rich structure of the scalar sector. By doubling the doublets 
or/and adding the gauge singlets such next-minimal extensions of SM model have been con- 
sidered and studied from various perspectives. Clear motivations for such extensions is to 
reduce some SM shortcomings , like the better agreement with precision electroweak fit, the 
theoretical problem of mass hierarchy or the dark matter problem. 

Having more interacting scalar bosons one can expect qualitative changes in the scalar 
boson sector. In some circumstances the binding forces between scalars can appear strong 
enough to produce the bound states. What is the spectrum of appearing bound states and 
how they exhibit their existence in the collider experiments are important questions. Most 
actually , with the run of LHC, what do we observe if the Higgsonium realizes in the nature? 

As a simplest model choosen for the actual calculation is the extension of SM that in- 
volves the addition of a real scalar singlet S to the SM Lagrangian. The phenomenological 
implications for singlet extension SM (xSM) has been studied from the collider and cosmo- 
logical perspectives [l|. The later typically require small mixing with the SM Higgs and 
from the perspective of bound states it reduces to the SM. In such circumstances it was 
shown in [2j that super heavy Higgs wlh — ITeV would be needed to form a bound state. 
In our model we will consider large mixing which could produce two scalars Hi and H 2 both 
having masses at few hundred GeV. 

II. MODEL 

In what follows I will consider the xSM, the SM is obtained by putting the new couplings 
to zero. The Lagrange density for the xSM model is 



where H denotes the complex Higgs doublet and S the real scalar. The term linear in S is 
chosen to vanishes after the spontaneous breaking. The potential is given by 



In unitary gauge the charged component of the Higgs doublet H becomes the longitudinal 
components of the charged H^-bosons and the imaginary part of the neutral component 
becomes the longitudinal component of the Z-boson. In Unitary gauge the Higgs field 
doublet reads 




(2.1) 




(2.2) 




(2.3) 
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Within our notation the mass terms in the scalar potential become 

Knass = \ i/ h h 2 + /4^ 2 + fiahS) , (2.4) 

where 

fi 2 h = 2\v 2 , /4 = K 2 + 8 2 v 2 , fi 2 hS = 25iv. (2.5) 

The mass eigenstate fields H\ >2 are linear combinations of the Higgs scalar field h and the 
singlet scalar field S. Explicitly, the inverse transformation reads 

h = c Hi — s H 2 ; 

S = sH x + cH 2 , (2.6) 
where c = cosO, s = sinO and the mixing is determined as 0] 

tanfl = j -, x = ^^-. (2.7) 

for positive x and for heavier singlet (x < 0) we have for mixing angle 



1 + vT 



x 
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tan9 = — . (2.8) 

\x\ 

The terms in the scalar potential that break the discrete S — > —S symmetry are propor- 
tional to the couplings 8\ and k% and we do not consider these term should very small we 
assume they are large enough to make a bound state and sufficiently strong communication 
with the rest of the Standard Model. The small mixing scenario with the light singlet like 
scalar has been considered in Q]. 

In this paper we will consider relatively large mixing angle 9 with Higgsiess masses less 
then 200 GeV which could be promising for experimental observation thorough LHC era. 
For such case the constraints from electroweak precision observables and their implications 
for LHC Higgs phenomenology has been already analyzed in [l|. 

The mass eigenstates satisfy in any case 

Ml = fi 2 h c 2 + fi 2 s s 2 + n 2 hS cs 

Ml = f j, 2 s 2 + l i 2 c 2 - l i 2 s cs. (2.9) 
The rest of the scalar potential contains four new parameters which are added to the SM. 

Vmt = + ^±s 4 + Xvh 3 + ^S 3 + 6 -^h 2 S 2 + 6 -±h 2 S + 6 -^hS 2 (2.10) 
4 4 3 2 2 v 

Based on the nonrelativistic consideration, the cubic interaction should be sufficiently en- 
hanced against the quartic one, otherwise the bound state cannot be formatted. 

The purely cubic interaction between mass eigenstates can be written in the following 

way 

Vcub = gmHl + gmH?H 2 + gmHiH* + g 222 H\ (2.11) 
with the new couplings g's are known cubic polynomial of sin# and cosO. 
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It is advantageous to evaluate new quartic couplings between physical states : 

V 4 = 9 -fHl + 9 -±fHlH 2 + 9 -±fHlHl + *f-H x Hl + # 2 4 (2.12) 

where A' are known quartic polynomial of cos and sin of mixing angel. 

Higgsonium as any other two body states in Quantum Field Theory is described by two 
body BSE: 

r = / VG [2] T (2.13) 



k 



where we use shorthand notation J k = i J and is the two particle propagator of the 
constituent Higgsies H\. In momentum space it can be conventionally written as 

G [2] {k,P) = D{k + P/2,M 2 )D(-k + P/2,M 2 ); (2.14) 

Let assume that the attractive interaction between heavy Higgsiess H 2 is strong enough 
to form a bound state. Within the xSM the irreducible BSE kernel in lowest order reads 



V = 6A im + E 



x=s,t,u 



40ii2 + 36<7? u 



x-M£ x- M{ 



(2.16) 



where the first term represent the pure constant interaction and the s, t, u are usual Mandel- 
stam variables. Down indices show which field - Hi or/and H 2 belong to a given interaction 
vertex. 

It is advantageous to explicitly divide the solution which is independent on the relative 
momenta of the constituents. Following the notation in 0] the original BSE can be rewritten 
in the following form 

r p ( P ,p) = r,(P) / v p (k, P ,p)G^(k,p)+ [ v p (k, P ,p)G®(k,p)r p (k,p), (2.17) 

where 

Vi = V c + V.,V p = V t + V u , (2.18) 

The first term is supposed to collect all constant term, i.e the ones the do not depend on 
the relative momenta. In our tree kernel approximation reads 



4ff 2 2 21 36g 2 2 
P 2 - Ml P 2 - Mf 



Vi = 6A im + -^h~ 2 + -^TH (2.19) 



where the full solution of BSE is given by the sum 

r(p,p) = r I (P) + r p ( P ,p) (2.20) 

The equation for the function Tj(P) is purely algebraic 

TliP) ~ l-Vjf.G^P) (2 - 21) 



The BSE represents the singular equation which can be solved by some known method. 
The one known possibility is to perform a Wick rotation for a relative momenta of con- 
stituents while keeping the total square of four momenta P 2 timelike. 

The other well known possibility is the Minkowski solution performed within the utiliza- 
tion of the unique integral representation of the kernels and amplitudes that appear in the 

bse 

The bound state vertex function can be expressed as the following 



n P, p) =U]to £££ (2.22) 



-1 



where n is a integer and all the singularity is trapped by the zeros of the denominator in 
flP2D which reads 

F{a,r)] P,p) = a - (p 2 + P.pz + — ) - ie (2.23) 

Recall the known property of superrenormalized model studied yet: The function is 
more smooth as larger n is used. In practise, the BSE were solved for generalized Wick- 
Cutkosky models for the lowest value only, i.e. n = 1,2. The studied models were very 
simple ones with the most cubic scalar interaction presented. 

Here, having the quartic interaction presented, the inhomogeneous term is generated 
represented by Tj(P) which is just real constant for a given discrete value of bound state 
mass P. To avoid more complicated distribution we naturally assume 



1 LXJ 

r(p, P ) = r/(p) + J d v J da 



[F(a,r);P,p)} 



(2.24) 



where p p (a,rj) is assumed to be a real function, and not delta distribution. It fully corre- 
sponds to the function T p , noting its structure is fully driven by pure triplet interaction of 
Higgs. Furthermore, we explicitly choose n — 1 in the IR f)2.24p . following the most easy inte- 
gral representation of the inhomogeneous term in the expression, i.e. f k V p (k,p, P)G^(k, P). 
These integral corresponds to the Feynman scalar triangle diagram. 

One can show that the BSE can be converted to the regular integral equation for p p . It 
reads 



1 



a - Mf 



1 CXJ 

Ti(P)pi(a, 7)) + J dz J da p p (a, z)V( 



a, rj, a, z) 



'2.25) 



where pi,V are known regular function and the constant term arises due to the quartic 
interaction 



r/(P) 



Vj J dz J dap p (a, z)I F (P ,a,z) 

— 1 — oo 

i - vp-if^P 2 ) 



(2.26) 
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where V R is the renormalized constant interaction and If, Ib stand for the known one loop 
triangle and one loop bubble integrals. 

To renormalize the momentum subtraction renormalization scheme with zero momentum 
scale is used, 

Vi - 6X ^ + pT^Ml + p^TM2 (2-27) 
where Af m is the renormalized quartic coupling of heavier Higgs mass eigenstates H±. 



III. RESULTS 

After a suitable normalization then the BSE for the weight function (I2.25P has been 
solved by the method of iteration. The coupling constants have been varied to reach the 
real discrete spectrum of Higgsonia. 

Firstly I mention the result for SM Higgs. The only known input is the Higgs vev. 
v = 275GeV while the Higgs mass and the cubic coupling depends on the experimentally 
unknown A, these satisfy = y/2Xv, A3 = 2Xv. There are no bound states bellow certain 
critical coupling of A3. The first Higgsonium of the mass M = 2m^ is formed when = 
1.3TeV. Such a result doesn't provide reliable answer, moreover it brings more questions. 

Such a heavy Higgs boson is ruled out by the electroweak precision test. Furthermore, for 
such a fat Higgs, the Higgs sector of SM represents strongly coupled field theory and our BSE 
becomes only rough estimate of reality. In addition, switching on the top quark Yukawa, 
the fat Higgs becomes broad resonance and its fast decays should prevent the formation for 
bound states. 

To have a reasonable model which is not completely ruled out by electroweak oblique 
correction constraints we assumed relatively light scalars in the xSM. There is no bound 
state unless the new cubic coupling is large enough. To compare various quantities, the all 
dimensionfull quantities are scaled in units of Higgs vev. v. Then roughly pronounced, the 
new cubic coupling must be several times larger then the rest of Lagrangian parameters. This 
is the main conclusion from the numerical inspection of the large regime of xSM parameter 
space. 

Here I present the first preliminary numerical result: The parameters used as an input 
were as the following: v = 275GeV, A = 0.20, 5 1 = 1.20 * v, 5 2 = 0.40, k 2 = 0.10 * v, k 3 = 
5.0* v, k a = 0.20. 

It leads to two massive eigenstates with Mi = 179. 5GeV, M 2 = 177.7GeV, and couplings 
gm ~ 280GeV, g 2 22 — 4:00GeV , Ann — 0.33 and the appropriate mixing cos8 = 0.696. 
Solution of BSE gives the bound states 20 % lighter than the Higgs production threshold : 

M B = 0.8 * 2 * Mi = 286GeV . (3.1) 
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